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Does a prestellar core always become protostellar? Tracing the 
evolution of cores from the prestellar to protostellar phase 
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ABSTRACT 

Recently, a subset of starless cores whose thermal Jeans mass is apparently overwhelmed by 
the mass of the core has been identified, e.g., the core L183. In literature, massive cores such 
as this one are often referred to as "super- Jeans cores". As starless cores are perhaps on the 
cusp of forming stars, a study of their dynamics will improve our understanding of the trans- 
ition from the prestellar to the protostellar phase. In the present work we use non-magnetic 
polytropes belonging originally to the family of the Isothermal sphere. For the purpose, per- 
turbations were applied to individual polytropes, first by replacing the isothermal gas with a 
gas that was cold near the centre of the polytrope and relatively warm in the outer regions, 
and second, through a slight compression of the polytrope by raising the external confining 
pressure. Using this latter configuration we identify thermodynamic conditions under which a 
core is likely to remain starless. In fact, we also argue that the attribute "super- Jeans" is sub- 
jective and that these cores do not formally violate the Jeans stability criterion. On the basis 
of our test results we suggest that gas temperature in a star- forming cloud is crucial towards 
the formation and evolution of a core. Simulations in this work were performed using the 
particle-based Smoothed Particle Hydrodynamics algorithm. However, to establish numerical 
convergence of the results we suggest similar tests with a grid-scheme, such as the Adaptive 
mesh refinement. 
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1 INTRODUCTION 

Stars form out of relatively dense volumes of molecular gas, typ- 
ically denser than ~ 10 4 cirT 3 , generally referred to as cores. The 
dynamical stability of a core is determined by the relative import- 
ance of thermal (which includes turbulence), and magnetic pressure 
against its own gravity. A "starless" core is one in which there is no 
evidence of star-formation. A gravitationally bound starless core 
with internal temperatures typically of order a few Kelvin, and a 
flat central density relative to the protostellar core, is called a pre- 
protostellar or simply prestellar core (Ward-Thompson et al. 1994); 
see also Di Francesco et al. (2007) for a relatively recent review of 
the observed characteristics of starless and prestellar cores. 

Having become sufficiently massive by accreting gas from the 
ambient medium, a core may evolve to become centrally dense, 
mimicking the density profile of the singular isothermal sphere. 
The protostellar phase of a core is thus one in which a subsequent 
collapse leads to the formation of a young stellar object, often de- 
tected due to its emission at infrared wavelengths from warm dust 
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surrounding the young protostar. Young stars also launch outflows 
and therefore, are useful tracers of star-formation activity. 

The process through which a core becomes protostellar is not 
fully understood, though a possible evolutionary track for prestellar 
cores was proposed by Simpson et al, (2011). This process, how- 
ever, may be crucial towards explaining the distribution of stellar 
masses characterised by the initial mass function (e.g., see Ward- 
Thompson et. al. 2007). In general, prestellar cores should collapse 
when internal support such as thermal or magnetic pressure is over- 
whelmed by gravity. Several cores that appear to violate this simple 
evolutionary model have been identified, however. For instance, the 
core L694-2, though starless, shows signs of infall (e.g., Harvey et 
al, 2003), many others show signs of outward motions or both in- 
ward and outward motions, like B68 (Lada et al. 2003). 

In addition, several cores have been identified as being star- 
less despite having either: (i) masses apparently greater than their 
thermal Jeans mass, for a non-magnetised core (e.g., see Sadavoy et 
al. 2010a,b), or (ii) relatively low internal velocity dispersion, of or- 
der ~0.1 km/s to 0.3 km/s with cold interiors (e.g., cores L1498 and 
L1517B; Tafalla et al. 2004). Cores of type (i) are also referred to as 
"super- Jeans" in literature. Including line observations to determ- 
ine velocity dispersion, Schnee et al. (2012) recently determined 
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that the continuum-derived masses of a sub-sample of five super- 
Jeans cores in Ophiuchus and Perseus are about equal to the virial 
mass of a critically stable Bonnor-Ebert sphere, suggesting weakly 
that the super- Jeans cores are bound even when non-thermal mo- 
tions are considered. The super- Jeans nature of these starless cores 
is of course debatable and will be examined in this contribution. 

Observed lifetimes (and estimated long lifetimes of the order 
of several free-fall times), of some prestellar cores suggest they 
are not in free-fall, and perhaps supported by some mechanism. 
Magnetic fields and turbulence may play a key role in this regard. 
For instance, cores located in magnetised regions of a cloud should 
evolve on the ambipolar diffusion timescale, usually of the order 
of a few Myrs, although cores close to super-criticality are likely 
to evolve much faster, on the order of only a Myr (e.g., Ciolek & 
Basu 2001). The age of a core can be estimated through a study of 
the gas-phase abundances of molecular species such as CS, NH 3 , 
CO and its isotopologues, and other complex molecules. Ages so 
derived for the cores L1498 and L1517B (Tafalla et al. 2004), of the 
order of 1-2 Myrs are indeed, consistent with the ambipolar times- 
cale for weakly magnetised cores. Similarly, the more chemically 
evolved cores in the Pipe Nebula lie in strongly magnetised parts of 
the cloud, and probably will also evolve on the classical ambipolar 
diffusion timescale, of order ten free fall times (Frau et al, 2010). 
According to numerical simulations, however, turbulence, the other 
possible support against gravity, decays on a relatively short times- 
cale, typically less than a sound-crossing time (e.g., Pavlovski et al. 
2002). Consequently, it is unlikely to explain the long lifetimes of 
some prestellar cores. These cores have therefore been, the sub- 
ject of intensive investigation by several authors, e.g., Evans et 
al. (2001), Crapsi et al. (2005a,b), Hatchell et al. (2007), Ward- 
Thompson et al. (2010), to name only a few. 

Some authors have studied the formation of cores numerically 
via the growth of density perturbations. Vazquez-Semadeni et al. 
(2005) reported that some of these perturbations condensed out as 
small clumps, a few of which collapsed while others re-expanded 
before eventually dispersing. Similar conclusions were also derived 
by Klessen et a/.(2005) from their simulations. Galvan-Madrid et 
al. (2007), simulating turbulent, magnetically super-critical clouds, 
derived results consistent with the findings of the latter authors. 
In fact, cores forming in the simulations by Galvan-Madrid et 
a/. (2007) evolved on a relatively short timescale, typically 1-2 
Myrs, as suggested for the cores L1498 and L1517B. The well- 
known starless core B68 on the other hand is probably in a critical 
equilibrium and might perform low-amplitude oscillations when 
perturbed weakly (Keto et al. 2006; Redman et al. 2006). In the 
present study, we will examine the stability of cores and determ- 
ine new conditions that are likely to retard them from collapse. The 
plan of the paper is as follows. We commence in §2 by briefly dis- 
cussing the problem, followed by a short recap of the physical prop- 
erties of the Isothermal sphere and an introduction to the numerical 
scheme used here. Simulations are discussed in §3 and §4 before 
concluding in §5. 



2 INITIAL CONDITIONS 
2.1 Description of the problem 

The possibility of the formation of a prestellar core out of a natal 
volume of gas will be examined in this work. This hypothesis will 
be examined by running a set of comparative models that will use 
two varieties of the Isothermal sphere. Note that, "Isothermal" with 



BONNOR-EBERT SPHERE 




COLD CENTRAL REGION 



Figure 1. A schematic illustrating the nomenclature adopted in the rest of 
this paper. Shown here is a core of radius, r care , assembled by radial inflow 
of gas. The core has been shaded grey and encloses the cold central region 
shaded black, while the outer annulus in a lighter shade represents the poly- 
trope, a Bonnor-Ebert sphere (a variety of the Isothermal sphere used in 
some simulations here). 



"I" will be reserved for the polytrope generated by the solution to 
the Lane-Emden equation, Eqn. (2) below; "i", will be used to de- 
note uniform gas-temperature. The isothermal gas within the re- 
spective polytropes will be superposed by a more realistic temperat- 
ure distribution, often found in prestellar cores and defined by Eqn. 
(7) below. The resulting assembly is demonstrated by the schematic 
sketched in Fig. 1, where r core is the radius of the core assembled at 
the centre of the agitated polytrope. A core in this work is defined 
as the pocket of gas within the model volume having average dens- 
ity greater than a threshold, p,/„ e ,/, = 10~ 19 g cirT 3 ~ 10 4 cm 4 . 
The polytrope is perturbed by inwardly propagating disturbance 
triggered due to a slight increase in the external pressure, achieved 
by suitably raising the temperature of the gas confining the poly- 
trope. 



2.2 The Isothermal sphere revisited 

The principal equation characterising a polytrope is the well-known 
Lane-Emden equation, 



1 da 



(i) 



Gas pressure, P, under the isothermal approximation, is related to 
its density, p, via a simple equation of state of the form, P = a\p y , 

where oq = j , is the sound-speed for isothermal gas held at 
temperature, To- The exponent, y = 1 + 1/n, where the polytropic 
index, n = oo, for an isothermal gas and the mean molecular mass, 
m = 4 x 10~ 24 g, for a typical cosmic mixture (i.e., a mixture of 
molecular hydrogen and helium). 

The Isothermal sphere belongs to a larger family of polytropes 
defined by Eqn. (1). Under the isothermal approximation Eqn. (1) 
reduces to 



df.1 



(2) 



where fi = |f 2 ^|j, the Lane-Emden function, if/ = , and £ 

is the dimensionless radius (Chandrasekar 1939). Stability charac- 
teristics of the Isothermal sphere show that it is inherently unstable 
against gravity for any radius, £ B = £ > = 6.45. The dens- 
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Table 1. List of models tested. 



Serial Confinement Polytropic Temperature 
No. 



Observation 



model 


distribution (T(r)) 








Pt 


Equation 7^, T in = 6 K, T out 


= 15 K 




Collapse 


A 


Equation 7, T m = 6 K, T imt = 


15 K, Tjcm = 


50 K 


Collapse 


B e 


Equation 7, r,„ = 6 K, T out = 


15 K, Ticm = 


20 K 


Collapse 


B 


Equation 7, r,„ = 10 K, T oul 


= 35 K, Ticm 


= 40K 


Oscillations 


B 


Equation 7, T\ n = 6 K, T out = 


35 K, Ticm = 


40 K 


Oscillations 




T = 20 K, T ICM = 25 K 






Collapse 



c 
c 
c 
c 



a U-Unconfined polytrope 

^ Originally a singular Isothermal sphere but now superposed with a radially dependent temperature profile, T(r) 

c T(r) defined by Equation (7) 

d C-confinement by external pressure defined by Eqn. (3) 

e Originally a Bonnor-Ebert sphere but now superposed with a radially dependent temperature profile, T(r) 

t Isothermal 



ity, p(^), of the Isothermal sphere decreases radially outward and at 
f = £crii, it falls to about 70% of its value at the centre. 

Configurations with £ ^ £.„■,, called a Bonnor-Ebert (BE) 
sphere, is confined by a finite external pressure, P ext , and supported 
internally by thermal pressure against gravitjQ. The pressure con- 
fining a BE sphere is 



Pext = Pee 



(3) 



where p c is the central density, and the asymptotic solution to 
the Isothermal Lane-Emden equation (Bonnor 1956, Ebert 1955). 
By demanding pressure equilibrium at the boundary of the Bonnor- 
Ebert sphere the minimum temperature, T , necessary to maintain 
stability can be readily obtained by combining Eqn. (3) with the 
equation of state for an isothermal gas. Thus, 



(4) 



k B is the Boltzmann constant. Finally, the maximum mass within a 
BE sphere of radius £ s is, 



M BE = M{% = fr) = 4npX^ B X 



(5) 



/ « 2 \ 1/2 

where the length scaling factor, R = I 1 , and p(^) has been 
defined above for Eqn. (2). 

2.3 Smoothed Particle Hydrodynamics (SPH) 

The well-known Lagrangian, particle-based scheme SPH finds a 
wide-range of applications in computational hydrodynamics, in- 
cluding astrophysical hydrodynamics (e.g., Monaghan 2005). SPH 
particles are not point masses in the strictest sense so that each 
particle has finite spatial extent defined by its smoothing length, hj, 
apart from other state variables, viz. mass, velocity and temperat- 
ure. Normal SPH particles exert the gravitational and the thermal 
force on each other and forces are calculated at every time step. 
After-ward, the position of each particle is revised by integrat- 
ing the equation of motion. SPH particles interact with each other 



I V' 2 

1 The free-fall time, Iff = , defined for gas with average density, 

p, though strictly inapplicable to a polytrope, has been adopted as the unit 
of time in the rest of this article for mere convenience. 



through a parametrised numerical viscosity that prevents them from 
penetrating each other. Viscous dissipation becomes effective for 
pairs of particles approaching each other and diminishes rapidly 
for those receding. 

Computational efficiency is acquired by stacking particles on 
an octal tree, like the one suggested by Barnes & Hut (1986), where 
leaf cells at the bottom of the tree contain less than four particles. 
Calculation of the net force and the search for nearest neighbours, 
N ne n, s , of a particle are carried out with this tree. Though contri- 
butions from distant cells are approximated, the accuracy of cal- 
culations is improved by including the quadrupole moments. The 
time integration is performed by using the multiple time-stepping 
algorithm where each particle is assigned an appropriate time step. 

The temperature of gas particles is calculated using the radi- 
ative transfer scheme introduced by Stamatellos et al. (2007). The 
scheme proceeds by constructing a polytrope typically of index, 
n=2, around each SPH particle, and then calculates the average op- 
tical depth for that particle. The optical depth determines the ex- 
tent to which a particle is shielded and calculates the correspond- 
ing change in internal energy of that particle. The algorithm uses 
parametrised opacity, calculated such that it accounts for physical 
processes that dominate a density-temperature range. See Table 1 
in Stamatellos et. al. (2007) for a list of processes covered. Fi- 
nally the revised temperature of that particle is calculated by solv- 
ing the equation for thermal equilibrium. Gas is assumed to be 
primarily composed of hydrogen(70%) and helium(30 %), even as 
the scheme successfully mimics the relevant gas-phase chemistry 
within a typical prestellar core. The numerical code SEREN used in 
this work is a well-tested algorithm; see Hubber et al. (2011) for 
further details. 

An SPH particle denser than a predefined threshold, 
Psink ~ 10~ 15 g cnT 3 , is replaced by a sink if it also satisfies two 
other qualifications : (a) it has the desired number of neighbours, 
at least N„ eihs , and (b) it is gravitationally bounded (Bate et al. 
1995). The sink particle in the present set of simulations represents 
effectively the protostellar phase of a putative star-forming core. 
The reason for choosing a relatively low p sj „ k will become apparent 
below. This choice of p si „t also implies that temperature of the 
gas within our test polytropes exhibits little change throughout 
a simulation as adiabatic heating is ineffective below ~ 10~ 12 g 
cm 4 . Simulations discussed here were performed on a 12-thread 
dual-Xeon X5675 processor of the Hydra super-computing cluster 
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Figure 2. Shown in the left- and right-panel is the respective radial density profile for a singular Isothermal sphere (fg=10) and a Bonnor-Ebert sphere(fs=4), 
each of unit mass and radius; the plot is made in dimensionless units. Marked by dots is the SPH calculation of density while result from analytic calculation 
is plotted by a continuous line. 



at the Indian Institute of Astrophysics. 

Resolution Simulations were performed with 5.5 x 10 6 particles, 
of which the number of gas particles, N gas = 5 x 10 6 , while the 
rest were particles representing the intercloud medium (ICM). 
The initial average SPH smoothing length, h mg , that provides an 
estimate of the spatial resolution is defined as 



h 



R l N \ 1/3 
"mtf / 1 'ueibs \ 



neibs \ 



(6) 



where A', 11 .,; M =50. Thus, h mg ~ 1.5 x 10~ 3 pc, so that, X = ro/h avg ~ 
15, i.e., the central region within radius, r , is resolved by about 10 
SPH particles; r is defined by Eqn. 7 below. Increasing X by an 
order of magnitude, would require N gas to be at least 3 orders of 
magnitude higher than the current choice, well beyond our present 
computational capabilities. The current choice of N gas , as will be 
seen below, is sufficient, for our extant purpose since we are not 
keen to resolve details of the protostellar phase. 

The minimum resolvable density, p ms ~ 1.1 x \Q~ 20 (N gas ) 2 
g citT 3 , which for the parameters chosen in the present work is 
~ 1CT 10 g cm 4 , about five orders of magnitude greater than p S i„t. 
Similarly, the minimum resolvable mass, M res , is 



M re 



2N„, 



where the bracketed quantity in the expression above is simply the 
mass of a single SPH particle, m, (Bate & Burkert 1997). The Jeans 
mass, Mj, at the sink-formation threshold of p sink ~ 10~ 15 g cnr 3 
and temperature 10 K (see also Table 1), is ~ 1.6 x 10~ 2 M G . Thus, 
■j^ — 1500, implying, a Jeans mass will be resolved by about 1500 
particles. Our simulations are therefore well resolved to study the 
onset of the protostellar phase, and can even be extended to the 
formation of the first core. However, the present study does not 
proceed beyond the protostellar phase. 



2.4 The test polytropes 

Out of the infinitely large family of the Isothermal sphere, we ini- 
tially opt for one each with radius on either side of its critical ra- 
dius, with £ s = 4 and f g = 10, respectively. The former is a 
Bonnor-Ebert sphere while the latter, a singular Isothermal sphere 
that is inherently unstable to gravity. The initial density profiles for 



the respective configurations are shown in Figure 2. Details about 
assembling the respective polytropes have been described in the 
following subsection, §2.5. 

Having assembled each of the two settled polytropes, a non- 
isothermal temperature distribution, T(r), of the type 



T(r) = T oul - 



{Tout T jjj) 

i + oT^oF 



(7) 



initially suggested by Crapsi et al. (2007) for the starless core L1544 
was overlaid. Here T mt = T(r > r ), T in = T(r < r ), and r ~ 
5000 AU. This expression accounts for the relatively low temper- 
atures within the interiors of prestellar cores, which according to 
detailed radiative transfer modeling, seldom exceed 10 K. For at 
typical prestellar densities the gas-dust coupling becomes domin- 
ant and dust particles cool efficiently via thermal radiation even as 
gas-phase coolants freeze out onto cold dust grains (e.g., Goldsmith 
2001). The external regions though, remain relatively warm due to 
a combination of a number of heating sources like cosmic rays, en- 
ergetic external photons, turbulence, etc. While this is the case for 
real cores, in the present work however, there is no external source 
of heating. 

The initial mass, M,„„, and radius, r,„„, of the test polytrope 
in the present work is respectively, 5.5 M Q and 0.15 pc. To main- 
tain stability for the above-defined choices of mass and radius, gas 
within the test polytropes will have to be maintained at To > 5 1 K, 
according to Eqn. (4) . However, our choice of the gas temperature 
renders the polytropes dynamically unstable; see Table 1. 



2.5 Assembling the polytropes 

Polytropes used in the actual simulations were generated by first as- 
sembling two Isothermal spheres of unit mass and radius, one each 
for £ s = 10 and 4, respectively. We remind, is the dimensionless 
radius of the Isothermal sphere. Polytropes of either type were as- 
sembled by positioning particles randomly within them after which 
they were settled, before finally scaling them to the desired dimen- 
sions. A settled polytrope is one in which the Poisson noise asso- 
ciated with the random assembly particles has been significantly 
dissipated. This was achieved by repeatedly revising particle posi- 
tions within the polytrope in the presence of artificial viscosity that 
dissipates thermal noise. The polytrope was prevented from diffus- 
ing into vacuum by jacketing it with a warm intercloud medium 
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Figure 3. The coeval plots in the left- and right-hand panel show respectively the radial profile of the gas-density, and the radial component of gas-velocity 
within the polytrope A for case 1. Radial distance is marked in units of the precollapse radius, r,,„ f , of the polytrope. Evidently, as the interiors of the polytrope 
collapse, the gas starting from rest falls-in increasingly rapidly although, it always remains subsonic (sound speed within the envelope of the precollapse 
polytrope is 0.23 km/s). 



0.06 t>c ► 



og density 

Figure 4. A rendered density plot showing the mid-plane close to the centre 
of the polytrope in case 1 at a latter stage(t ~ 0.3 Myr =1.1 tft). Overlaid on 
this plot is the local velocity field shown by arrows, while the sink particle 
representing a young protostar has been marked by an * at the centre of the 
singular polytrope. The radially inward moving gas within the polytrope 
can be readily inferred from the direction of the velocity vectors. 



(ICM), represented by a special type of SPH particle that only ex- 
erts hydrodynamic force on gas particles. The polytrope-ICM sys- 
tem was then confined by another thin shell of boundary particles, 
meant to prevent SPH particles from escaping into vacuum. These 
boundary particles remain fixed in space and do not exert any force 
on other SPH particles. 

The number density of particles was adjusted so that pressure 
equilibrium was maintained across the gas-ICM interface. Figure 2 



S. 0.03 




Figure 5. The plot shows the early accretion history of the young protostar 
in case 1, represented by a sink particle in the simulation. 



shows the density profiles for settled polytropes of unit dimensions. 
Having so settled the two polytropes, they were then rescaled to 
the desired dimensions, and a radially dependent temperature pro- 
file defined by Eqn. (7) above was then superposed on either of 
them. Thus, the polytropes used in actual simulations, with the ex- 
ception of case 6 listed in Table 1, were no longer isothermal. In 
the rest of this article the polytrope with radius £; B = 10, originally 
a singular Isothermal sphere, but now a singular non-Isothermal 
sphere, will be referred to as polytrope A. Similarly, the Bonnor- 
Ebert sphere with £ s = 4, also no longer Isothermal, will be referred 
to as polytrope B. Before proceeding to discuss the actual simula- 
tions it would be useful to demonstrate the suitability of the SPH 
algorithm for the extant purpose. The fact that the algorithm can 
reproduce the analytically predicted behaviour of the Isothermal 
sphere of critical radius, % B = £„-,,, is demonstrated in appendix A 
below. 



3 NUMERICAL SIMULATIONS 

We begin by introducing the convention adopted for velocity plots 
in the rest of this paper. Shown in these plots is the radial velocity 
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Figure 6. As in Fig. 3, radial profiles of the gas-density and radial velocity within the polytrope for case 2 have been shown respectively in the left-, and 
right-hand panels. While the plots are generally identical to those in Fig. 3 for case 1, gas close to the centre becomes mildly supersonic as indicated by the 
blue velocity curve (r ~ 1.1 ttt). 



of gas within a polytrope at different epochs. Gas moving towards 
the centre of a polytrope has negative velocity while that moving 
away from the centre has positive velocity. The direction of fluid 
motion within a polytrope is also reflected by the direction of ve- 
locity vectors overlaid on top of rendered density plots whenever 
necessary. 



over a free-fall time, slightly longer than that observed in case 1. 
The density and velocity plots for this case have been shown in Fig. 
6. As in the previous case, gas close to the centre of the polytrope 
in this case also starts from rest and moves inward, though after 
a free-fall time the gas-velocity becomes slightly supersonic (blue 
curve, t ~ l.lfyy); see Fig. 6. 



3.1 Cases 1 and 2 with polytrope A 

The singular polytrope A in cases 1 and 2 when perturbed, 
promptly collapses, a result that is hardly surprising. The polytrope 
for case 1 is the same as one in vacuum at infinitely large tem- 
perature. Thus, while gas within the interiors of this polytrope ( 
r l r init ^ 0.6 ~ 0.1 pc), is soon overwhelmed by gravity and begins 
to collapse, that in the outer regions diffuses. Starting from rest, the 
collapsing gas acquires a velocity ~ 0.2 km/s. Shown respectively 
in the left-, and right-hand panels of Fig. 3 is the gas-density and 
the velocity within the collapsing polytrope at different epochs for 
this case. 

Collapse of the polytrope effectively marks the beginning of 
the protostellar phase. Infall is also evident from the overlaid arrows 
that represent the direction of the local velocity field. The rendered 
density plot in Fig. 4 shows the mid-plane of the singular polytrope 
for case 1 where the young protostellar object that forms at the 
centre of the polytrope is marked by an *. Gas is steadily accreted 
by the young protostar and its mass, M„, grew almost linearly with 
time at a rate of a few times 10~ 5 M Q yr _1 , as can be seen from 
its early accretion history shown in Fig. 5. A comparison with the 
analytic estimate of protostellar accretion-rate, M„ would be use- 
ful, and the simplest approximation to which is the Bondi-accretion 
model. If p s and a are respectively, the average density of the gas 
being accreted and sound-speed in this gas then, M, = Anp s r]a\ 
r s - GI ^i" , is the sonic radius. Plugging in the appropriate values 
yields, M, ~ 3 x 10~ 5 M G yr" 1 , which agrees with that deduced 
from the early accretion history. This rate of accretion is also con- 
sistent with those reported for typical protostellar cores . Andre at 
al. (2007), for instance, have derived similar values for cores in the 
Ophiuchus star-forming cloud. 

The same polytrope as in case 1 , but now confined by a finite 
external pressure defined by Equation (3), was used in case 2. By 
contrast, the polytrope in this case is in global collapse that occurs 



3.2 Remaining cases with polytrope B 

We will now discuss the next set of our simulations, listed as 3, 4 
and 5 in Table 1 . Reminding that polytrope B is a pressure-confined 
BE sphere superposed with a non-isothermal temperature distribu- 
tion defined by Eqn. (7) would be useful. Results from this set of 
simulations are mutually contrasting, so while the polytrope in case 
3 became singular, it exhibited large-scale radial oscillations un- 
der different initial conditions in cases 4 and 5. Inferences drawn 
from this set of simulations will then be used in §4 to hypothesise a 
scenario in which prestellar cores could be arrested in their evolu- 
tion to the protostellar phase. As in the first two cases, gas initially 
at rest within the polytrope in case 3 gradually moves towards the 
centre of the polytrope as it is swept by an inwardly propagating 
compression-wave. The wave though, remains largely sub-sonic 
even as infalling gas close to the centre gathers momentum. This 
can be readily seen from the velocity plots on the right-hand panel 
of Fig. 7. 

The dynamically unstable polytrope in cases 4 and 5 failed 
to become singular and instead, performed large-scale radial os- 
cillations. As in the previous cases, the inwardly propagating 
compression-wave was again triggered through a slight increase in 
the confining pressure, P ext , Consequently, the polytrope showed 
intermittent phases of enhanced central density that eventually ex- 
ceeded the threshold, p,i, res i,, for a prestellar core. Rendered density 
plots overlaid with the local velocity vectors in the upper-panel of 
Fig. 8 show the compression-phase of the agitated polytrope, while 
the rarefaction is shown in the lower panel. These plots spanning a 
little over a free-fall time of the original polytrope readily demon- 
strate the assembly of a prestellar core due to the radial inflow of 
gas, assisted by the compression-wave (upper-panel, central-plate), 
followed by the first rebounce of the core (lower-panel, right-hand 
plate). Although the simulation was allowed to continue over six 
free-fall times, for want of space only a sequence over the first free- 
fall time has been shown in Fig. 8 here. 
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Figure 7. As in Figure 6, radial profiles of the gas-density and the radial component of gas-velocity within the polytrope for case 3 have been shown in the 
left- and right-hand panels, respectively. Polytrope B in this case is dynamically unstable and gas close to its centre, starting from rest, is soon overwhelmed 
by the inwardly propagating, subsonic compression-wave as can be seen from the velocity plots (sound-speed within the polytrope is ~0.23 km/s). 
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Figure 8. Upper-panel: A cross-sectional rendered density montage (in units g cirT 3 ) showing the first half-cycle of oscillation, where gas assisted by the 
compression-wave moves inward. As in Fig. 4, the local velocity field is indicated by arrows overlaid on the plots; t = 0. 1?//, 0.5f//, 1?//, respectively for plots 
in left-, central-, and right-hand panels (If// ~ 0.4 Myr). Lower-panel: The compression soon makes way for rarefaction, evident from the direction of velocity 
vectors in the plots on the left- (; = 1?//) and the right-hand panel (f = 1,2?//), respectively. These plots demonstrate the rebounce of the pressure-supported 
core. 



The rendered plots in Fig. 8 are further supplemented by those 
in Fig. 9 where radial variation of the gas-density and velocity 
within the agitated polytrope has been plotted for one oscillation, 
on the upper and central-panel, respectively. The gas-density within 
the polytrope steadily rises as the compression-wave moves in- 
ward, then having reached a maximum and as a rarefaction sets 
in, the density decreases (t ~ 1.2f//). Correspondingly, gas within 
the polytrope over-run by the inwardly moving compression flows 
towards the centre with ever-increasing velocity, as demonstrated 
by the curves corresponding to t = 0.5tff and / = 0.75?//. Then 
as a rarefaction sets-in, gas within the core assembled during the 



compression-phase begins to move outward gradually and its velo- 
city becomes less negative (/ = 1?//). The outwardly directed gas 
motion associated with the core-rebound is more explicitly visible 
from the velocity-curve for / = 1.2f//. 

Plotted in the lowest panel of Fig. 9 is the mass of the core, 
M core , and its radius, r core . The core assembled during the first half 
cycle of contraction and dissolved in the following rarefaction, is 
reassembled during the next half-cycle of contraction followed by 
rarefaction again. This sequence over six free-fall times of the ori- 
ginal polytrope is demonstrated by this plot. The radius of this core 
varies typically between 0.01 pc and 0.07 pc, comparable to core- 
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Figure 9. Shown on a logarithmic-scale in the upper-panel is the radial 
distribution of gas density within the BES for case 4 at different epochs. A 
coeval plot showing the radial component of velocity has been shown on 
the central-panel. The inward flow of gas assisted by the compression-wave 
is reflected by the negative gas-velocity. The compression-wave, as in the 
earlier three cases, remains largely sub-sonic, though gas close to the centre 
briefly becomes transonic (curve for t ~ 0.75f/y). Thereafter the core soon 
rebounds and gas moves outward as demonstrated by the characteristics for 
gas density and radial-velocity at t = lf//&1.2t//; see also text. The plots 
on the lower panel show the radius (red curve), and mass (green curve) of 
the assembled core. 



radii reported in several surveys of dense cores (e.g., see Jijina et 
al. 1999, Kirk et al. 2007, Belloche et al. 2011). In dimension- 
less units the core radius is, = ^g 5 -, where Rq is the physical 
scale-length introduced previously in §2.2. Interestingly, ^ re , for 
the assembled core, as defined above, is 2, confirming its pressure- 
supported nature. Its rebounce is therefore hardly surprising. That 
particles exceeding the density threshold, p,/„-„/, ~ 10~ 19 g cnr 3 , 
formulate the core would serve as a useful reminder for our readers 
at this point. 

The critical Bonnor-Ebert mass for the core, M c gi e , is obtained 
by simply rewriting Eqn. (5), now with = 2. This yields, 

M%r = M BE (^< = 2) ~ O-W^-L (8) 

The central density, p c , and temperature of gas within the core 
is ~ 10~ 17 g cnr 3 and 10 K, respectively so that M££ e ~ 0.032 
M Q . Now the mass enclosed within the radius, r cme = ^ re Ro ~ 
0.004 pc, is simply, M core (^ re ) ~ 4np c r 3 cor J3 ~ 0.04 M . Thus, 

M core (£ care ) 

J" = — — , is slightly in excess of unity. And further, it is 
significantly higher than unity when calculated using the mass 
within the entire core of radius, ~ 0.01 pc, pointing to its pro- 
clivity to collapse under gravity. Thus the plot in the lower-panel 
of Fig. 9 demonstrates the epochs when exceeded unity, i.e., 
portions of the green-curve exceeding M c °™ calculated by Eqn. (8) 
above. Crucially though, the core did not become singular. An im- 
portant caveat in the above argument is the assumption of a uniform 
temperature in Eqn. (8) to calculate the critical BE mass, M c °™ , for 
the core. This is inappropriate in view of the non-isothermal nature 
of gas within the test polytrope. In the strictest sense therefore, the 
stability criterion must be tested at different radii within the poly- 
trope. This is done by calculating the Jeans mass, Mj(r), and the 
BE-mass, M BE (r), as a function of the core radius. 

Thus, shown in Fig. 10 is the mass, M(r), within a shell 
of radius, r,(red characteristics), and the thermal Jeans mass, 
M,/(r2, (green characteristics), for these respective shells at differ- 
ent epochs. We note that the Jeans mass, Mjir), for a shell is pro- 
gressively lowered as the polytrope contracts although, it is never 
overwhelmed by mass within that shell, M(r). In other words, the 
compressional wave is not strong enough to assemble sufficient 
mass in a shell so that the Jeans stability criterion is never viol- 
ated for that shell. The sound-speed in the calculation of Mj(r) was 
replaced by the effective sound-speed to include the velocity dis- 
persion induced by the inwardly propagating disturbance. 

This plot is also supplemented by the one on the left-hand 
panel of Fig. B.l that shows the BE mass, M c °™((; = r/R ), as 
a function of radius after about 1 free-fall time just before the 
assembled core rebounded. Like the Jeans mass, Mj(r), the BE 
mass, Mg™ e (<f), for a shell also remains higher than the corres- 
ponding mass of the shell. Consequently, the polytrope does not 
become self-gravitating. This plot in conjunction with those in Fig. 
10 demonstrate that the Jeans criterion is essentially preserved at 
every radius within the polytrope, though Eqn. (8) leads to the con- 
verse conclusion. We may also add that the original result about the 
Super- Jeans nature of some starless cores by Sadavoy et al. (2010), 
was also the consequence of assuming an isothermal gas. So al- 
though the Jeans stability criterion appears to be violated, in the 
strictest sense it is not, as we have demonstrated here. 

To see if further reduction of gas temperature close to the 

2 Mj(r) = (ff/G) 1 5 (a 3 // (r)/p 1/2 (r)); a ef f(r), being the effective sound- 
speed for a shell of radius r, and average density, p(r). 
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Figure 10. Shown here plots comparing the mass, M(r), within a shell of radius r (shown in red), and the Jeans mass, Mj(r), for the respective shell (shown 
in green). Note that the plots are coeval with respect to the density and velocity characteristics shown previously in Fig. 9. 



centre of the polytrope led to a collapse, we expanded on case 4 
in the next case, listed 5 in Table 1 by lowering the central temper- 
ature to 6 K while maintaining a warm annulus of gas as in case 
4. The polytrope in this case also evolved in a manner similar to 
the previous case and the average temperature of gas within the as- 
sembled core was now 6 K. A plot of its gas-density as a function of 
radius of the polytrope at different epochs over one free-fall time is 
shown on the upper panel of Fig. 1 1 . Once again, it can be seen that 
having assembled a core at the centre of the polytrope, a rarefaction 
sets in and the central density is lowered. The velocity plots on the 
lower panel of this figure are therefore similar to those for case 4, 
plotted on the central-panel of Fig. 9. 

In the final case, listed 6 in Table 1, the polytrope B was made 
strictly isothermal at 20 K, well below the critical temperature, 
T , defined by Eqn. (4). As in all the other cases discussed before 
this, the polytrope was perturbed by slightly increasing the external 
pressure, P e „, that set off a compression-wave causing the initially 
gravitationally sub-critical polytrope to collapse and become sin- 
gular. Evidently, thermal support in this case was insufficient to ar- 
rest the collapse. Also as demonstrated by the plot in the right-hand 
panel of Fig. B. 1 , the critical BE mass, M^i^), is overwhelmed by 
the mass of the gas carried in by the compressional wave and so the 
polytrope becomes singular. Finally, shown in the left- and right- 
hand panel of Fig. 12 is respectively the radial density profile, and 
the radial component of gas-velocity within this polytrope. These 
plots are similar to those seen earlier for test cases 1, 2 and 3 where 



the respective polytropes also became singular; see also Figs. 3, 6, 
7, 9 and 11. 



4 DISCUSSION 

Using the two polytropic models A and B described in §2.4 above, 
we have attempted to reconcile some features of prestellar evolu- 
tion. Results presented in §3 demonstrate the importance of initial 
conditions towards the formation and further evolution of a pu- 
tative star-forming core and therefore determine its proclivity to 
spawn stars. Results from cases 1 and 2 demonstrate that the in- 
herently unstable polytrope A, aided by an inwardly propagating 
compression-wave, was easily driven to collapse. The polytrope B 
used in the remaining 4 simulations on the other hand, presents a 
very interesting case when perturbed similarly. The polytrope in 
cases 3 and 6 became singular as temperature of the gas within it 
was significantly lowered. 

On the contrary in test cases 4 and 5, where gas within the 
polytrope was maintained closer to the critical temperature, To, it 
exhibited intermittent phases of compression and expansion, a res- 
ult that we believe, could help improve our understanding of the 
relatively long-lived prestellar cores. While the compression- wave 
within the polytrope in these latter cases largely remained subsonic, 
the radial inflow of gas assisted by this wave assembled a pressure 
supported core at the centre of the polytrope in each of these two 
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cases. The core was perhaps only weakly bound. Having reached 
the prestellar phase (n ~ 10 6 cirT 3 ), this core rebounded, as can be 
seen from the radial density and velocity plots shown in Figs. 9 and 
11. 

Starless cores sometimes indeed, exhibit signs of inward mo- 
tion as also reported by Myers et al. (2000) for instance, and other 
references therein, for cores in the Perseus and Serpens molecu- 
lar clouds. Other similar examples include the well-known starless 
core L1544 (e.g. Tafalla et al. 1998), and a sample of those repor- 
ted by Lee, Myers & Tafalla (1999, 2001), Lee, Myers & Plume 
(2004), and more recently by Schnee et al. (2012). Inward-pointing 
velocity vectors within a core are generally considered to be pos- 
itive evidence for gravitational boundedness. However, in view of 
the results from cases 4 and 5, we suggest that while signatures of 
inward motion may be necessary to establish boundedness, they are 
certainly not sufficient. 

Having imposed a more realistic temperature distribution 
defined by Eqn. (7) on gas within the model polytropes in the first 
five of our test simulations, they were rendered non-isothermal. 
The polytrope for case 6 though isothermal, was cold at a temper- 
ature lower than a third of the that required to maintain dynamical 
stability. The dynamical evolution of polytropes observed here 
shows that the temperature within their interiors, i.e., within 
the radius ro, hardly affects their dynamical stability. For the 
compression-wave is set up in the outer regions of the polytrope 
and its strength depends on the temperature, T m „, of the gas in the 
annular region between the ICM and ro. We will revisit this point 
in §4. 1 where the possible implications for core-formation will be 
discussed. 

Other related work We emphasise that the compression-wave 
observed within either polytropic model, A and B, is predicted by 
the analytically derived self-similar solutions for the Isothermal 
sphere. A full set of self-similar solutions derived by Whitworth 
& Summers (1985), hereafter WS, show that those obtained 
earlier by Larson (1969), Penston (1969), Shu (1977), and Hunter 
(1977) belong to a larger family of solutions that can be broadly 
divided into 3 bands, with those respectively in bands and 1 
being natural solutions. These latter solutions by WS showed that 
a compression-wave over-running a SIS or a BES could possibly 
trigger a runaway collapse in their respective configurations, as 
is indeed observed here in cases 1, 2, 3, and 6. The collapse of 
the polytrope observed in these latter test cases agrees with the 
classic outside-in scenario of the Larson-Penston solution, and 
unlike the inside-out Shu solution. Similarly, solutions by Tohline 
et al. (1987) demonstrate that particularly strong fluctuations in the 
confining pressure could generate a compression-wave inducing 
the Isothermal sphere to collapse. 

The more interesting results from cases 4 and 5 though, fall 
somewhere in-between the natural-solution bands of WS. Either 
the inward-moving wave in these two cases is not strong enough, 
or conversely, the initial polytrope, though unstable, is sufficiently 
warm so that this wave remains largely subsonic as can be seen in 
the velocity plots of Fig. 9. The formation of a pressure-supported 
core prevents further collapse of the interior regions. The plots in 
Figs. 9 and 1 1 offer an instructive comparison with those in Figs. 
3, 6, 7, and 12. The agitated polytrope in these cases (4 and 5) is 
punctuated with phases in which the core mass appears to exceed its 
critical Bonnor-Ebert mass, M^ e , (see the plot in the lower panel of 
Fig. 9), and yet, it does not become singular. Although plots in Fig. 
10 contradict this conclusion. These latter plots demonstrate that 
the inwardly propagating wave is not strong enough to sufficiently 
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Figure 11. As in the upper-panel of Fig. 9, but for case 5. Shown here is the 
time- variation of the radial profile for gas density within the polytrope on a 
logarithmic scale is shown. 



raise the gas density, and therefore suitably lower the leans mass 
within the agitated polytrope so it could possibly be overwhelmed 
by the local gas mass. Consequently, there is no further collapse. 

This plot is further supported by another similar plot on the 
left-hand panel of Fig. B.l which shows that like the leans mass, 
the BE mass in case 4 also remains significantly higher than the 
gas mass at that point so that formally there is no violation of the 
stability criterion. It is thus evident that cores could possibly be 
classified as "super-Jeans" under a set of simplifying assumptions 
like isothermality of gas. In reality though, they are far from iso- 
thermal and follow a temperature distribution similar to the one 
defined by Eqn. (7) above. The commonly observed non-thermal 
velocity fields within starless cores can also be reconciled by the 
results from cases 4 and 5. 

Polytropic models have also been tested by other authors. 
Hennebelle et al. (2003), for instance, demonstrated that a Bonnor- 
Ebert sphere could be induced into collapse by sufficiently increas- 
ing the external pressure, P ext . Gomez et al. (2007), more recently 
discussed one-dimensional models where core-rebound was ob- 
served in some cases. Importantly, however, the rebound in their 
simulations was a result of the gravitationally sub-critical nature 
of the core-envelope system. In other words, the criterion for grav- 
itational stability was not violated in their simulations. Other au- 
thors, for instance, Ballesteros-Paredes et. al. (2003) have argued 
that cores forming in turbulent clouds are far from hydrostatic equi- 
librium and internal velocity dispersion is likely to delay their col- 



© 2002 RAS, MNRAS 000.[T]fT3] 



Does a prestellar core always become protostellar ? 11 



lapse. Turbulence though, is dissipative and loses energy rapidly on 
a timescale of the order of a free-fall time. It is therefore unlikely to 
support a core against gravity for too long. Furthermore, numerical 
simulations of turbulent cores show that, while some of the larger 
cores seem to produce filamentary structures, less massive cores 
produce protostellar accretion discs (e.g., Walch etal. 2010). Thus, 
while turbulence may indeed delay a collapse, it may not succeed in 
arresting one. Also, there is little any corroborative evidence since 
most starless cores show subsonic, or at best, transonic velocity 
fields as has been discussed above. 

It is clear that polytropic models are considerably successful 
in realising some dynamical properties of starless cores, though 
commenting on the chemical properties of a typical starless core 
is difficult. CO and its isotopologues freeze out at high densities, 
i.e., at densities in excess of 10 5 cirr 3 and temperatures less than 
~20 K (Di Francesco et al. 2007). Depletion of these species would 
be expected in centrally condensed cores and such depletion has 
been observed in the starless cores L1498 and L1517B (Tafalla 
et al. 2006). Although on the one hand there is depletion of CO 
species, on the other, there is enhancement of deuterated species 
as CO depletion curtails the formation of DCO + via the coupling 
between H2D + and CO (Crapsi et al. 2005b). Importantly, the core 
assembled in test cases 4 and 5 though ephemeral in nature, is suit- 
ably dense (n > 10 6 cm -3 ), and in fact, can perhaps sustain the 
associated gas-phase chemistry discussed above. 

Since there is not much variation of temperature during the 
expansion phase of the prestellar core, it is possible to reconcile 
apparent anomalies such as outward gas motions associated with 
chemical signatures of a typical dense starless core, or the con- 
verse. The starless core L429, for instance, has shown more red than 
blue asymmetry in observations of CS(2-1), indicating an outward 
motion (Lee et al. 2001). Line profiles indicative of simultaneous 
infall and outward motions, i.e., radial oscillations, across the face 
of the well-known starless core Barnard 68 were noted by Lada et 
al. (2003), and numerically modeled by Keto et al. (2006). How- 
ever, oscillations demonstrated by these latter authors had a much 
smaller amplitude in comparison to those demonstrated here in test 
cases 4 and 5. Furthermore, Lee & Myers (2011) in their recent 
study of starless cores suggested that early in their life, these cores 
might remain static, then perhaps enter an oscillatory phase before 
finally contracting, having become sufficiently massive. 

4.1 Implications for the formation and evolution of prestellar 
cores 

Supported by observations from two of our test models, 4 and 5 , we 
infer that radial inflow of gas can readily assemble a thermally sup- 
ported core. Such flows are more likely to occur in regions domin- 
ated by hydrodynamic instabilities such as the Kelvin-Helmholtz, 
Rayleigh-Taylor, or the Thin-shell instability which readily gen- 
erate filamentary clouds (e.g., Klessen et al. 2005; Anathpindika 
2009 a,b). See also Gong & Ostriker (2009), where the forma- 
tion of a core due to a supersonic, radial flow of gas was demon- 
strated. Furthermore, a core assembled in this fashion, depending 
on the strength of compression and gas-temperature, could either 
become singular or remain dynamically unstable over several free- 
fall times. In view of this latter possibility, we suggest that the star- 
lessness of some cores could be just a brief phase in the evolution- 
ary cycle of these respective cores. 

The contrasting nature of results from the set of simulations 
discussed here show that temperature of the gas within the model 
polytrope is crucial towards its dynamical stability. Results from 



the present work as well as arguments presented by Tohline et 
a/.(1987) re-emphasise the importance of cooling in potential star- 
forming gas. In star-forming clouds a differential temperature dis- 
tribution such as the one reported by Crapsi et al. (2005), and ex- 
ploited here could possibly arise due to variations in the efficiency 
of molecular cooling processes and local environment. Chemical 
processes and the environment responsible for cooling the gas dur- 
ing core formation therefore hold the key to the stability of that 
core. 



5 CONCLUSIONS 

Having accreted sufficient mass from its natal volume of gas a pres- 
tellar core becomes singular, spawning a young star within the col- 
lapsing core. Those cores in which this collapse appears to have 
either halted or seemingly reversed into expansion despite a cent- 
rally peaked density distribution may remain starless. In the present 
study we have attempted to explain these features with the help 
of a polytropic model that utilises two varieties of the Isothermal 
sphere, viz., the singular and the pressure-supported type. Both 
polytropes were allowed to evolve under self-gravity after super- 
posing a radially dependant temperature distribution. The eventual 
fate of either polytropes appears to depend critically on the initial 
distribution of gas temperature within it. 

A polytrope with cold gas is more likely to become singu- 
lar when perturbed by a disturbance of suitable strength. A similar 
perturbation, on the other hand, can also set up oscillations within 
the polytrope when gas within it is relatively warm. In the even- 
tuality such as the later, each cycle of oscillation was observed to 
last for a few times 10 5 years. The oscillating polytrope acquired a 
centrally peaked density profile even as density in the outer regions 
followed a relatively shallow power-law of the type, p(r) oc r~ k , 
similar to that observed in typical starless cores. We also report 
that in each of these contraction phases, mass within the central 
dense region, referred to as a core in this study, appears to exceed 
its critical mass for dynamical stability though it does not form- 
ally violate the stability criterion. Supported by this conclusion we 
suggest that the attribute "Super- Jeans" is subjective and could be 
potentially mis-leading. The assembled core survived for a few 10 4 
yrs before rebouncing during a rarefaction, only to be reassembled 
in the following cycle of compression. We therefore suggest, star- 
lessness is probably only a short, temporary phase in the lifespan 
of a pressure-supported core. Our simulations also offer a possible 
explanation for the origin of a non-thermal velocity field within a 
starless core and show that an inwardly directed velocity field is a 
necessary but not sufficient condition to establish its boundedness. 

However, it would be useful to repeat the present set of numer- 
ical tests with an Eulerian grid-based algorithm to further confirm 
the conclusions arrived at in this work. 

Without invoking the magnetic field or an external turbulent 
velocity field, the models discussed here reconcile the relatively 
long lifetimes of some cores. We have also demonstrated that pres- 
tellar cores can be assembled by a radial inflow of gas in star- 
forming clouds which is particularly likely in filamentary clouds. 
We also suggest that while inward pointing velocity vectors neces- 
sarily indicate inward gas motion, it is insufficient to establish the 
gravitational boundedness of a core with any degree of certainty. 
Extensive mapping of velocity fields within cores is essential to es- 
tablish if indeed their directions can be used as unambiguous pro- 
gnoses of the prestellar phase. It must be warned that although a 
polytropic model such as the one discussed in the present study 
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Figure 12. Coeval plots showing the radial profile of the gas-density and the radial component of the gas-velocity within the polytrope for case 6 have been 
plotted on respectively the left-, and right-hand panel. Either plots are similar to those seen previously for cases 1, 2 and 3 in Figs. 3, 6 and 7. The velocity plots 
show that as the compression-wave propagates inward, gas close to the centre, initially at rest, falls-in with ever increasing velocity as the polytrope becomes 
singular. 



appears commensurate with the dynamical properties of starless 
cores, it is not completely clear how in star-forming clouds such 
an assembly would be realised. The study of core formation and 
the effects of cooling processes on the local environment therefore 
deserve further attention. 
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Figure Al. Radial density profile of the Isothermal sphere with ifg ~ 6.45. 
The characteristics have been plotted at t ~ (black), and t ~ 0.3f// (red); 
tff ~ 1 Myr is the free-fall time of the Isothermal sphere. 

APPENDIX B: CAN A STARLESS CORE BE 
SUPER- JEANS ? 

The answer to this question is unlikely to be affirmative though 
some cores, such as the one assembled in this work,may appear 
to remain starless despite exceeding their thermal Jeans and/or the 
Bonnor-Ebert mass; they do not become singular. Plots in Figs. 10 
and B.l demonstrate that the compressional wave in the original 
polytrope in case 4 fails to assemble enough gas so that neither the 
Jeans criterion, nor the BE criterion for stability is violated and the 
respective polytrope fails to become singular. Instead it performed 
large amplitude radial oscillations and the assembled core reboun- 
ded. Under the assumption of isothermality though, the core ap- 
pears to have violated the stability criteria as shown in §3.2. This 
conclusion, however, is contradicted by a little more careful ana- 
lysis. 

This paper has been typeset from a TpX/ ETpX file prepared by the 
author. 



APPENDIX A: THE CRITICALLY STABLE 
ISOTHERMAL SPHERE 

This simple test was performed to demonstrate that the SPH al- 
gorithm could in fact reproduce the analytically expected singular 
nature of the Isothermal sphere having radius, % B = <f„-„ = 6.45. 
As in the main body of the paper, we first assemble and settle the 
Isothermal sphere, now for f s ~ 6.45. The settled polytrope was 
then rescaled so that it had a radius, ;„„, = 0.1 pc, and mass, M,,„, 
= 1 M Q . Gas within the polytrope was held at a uniform temperat- 
ure of 5 K (lower than the critical temperature, T , defined by Eqn. 
(4), which in this case is 14 K), surrounded by a jacket of particles 
mimicking the intercloud medium, maintained at 10 K. Without 
much surprise, the polytrope begins to collapse as can be seen from 
the radial density distribution plotted in Fig. Al at two epochs. 
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Figure Bl. Shown in these plots is the Bonnor-Ebert(BE) mass (green characteristic) as a function of the radial distance within the polytrope. Also shown is 
the corresponding mass, M(r)(red characteristic), at that radius. The plots in the left and the right-hand panel correspond respectively, to cases 4 and 6 after a 
free-fall time(~ 0.3 Myr). The core mass in case 4 fails to overcome its BE mass so that the polytrope fails to collapse. On the contrary, the polytrope in case 
6 becomes singular as mass close to the centre of the polytrope exceeds the corresponding BE mass. 
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